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Abstract 

We revisit the expansion recently proposed by Pulvirenti and Tsagkarogiannis for a sys- 
tem of N continuous particles in the canonical ensemble. Under the sole assumption that the 
particles interact via a tempered and stable pair potential and are subjected to the usual free 
boundary conditions, we show the analyticity of the Helmholtz free energy at low densities 
and, using the Penrose tree graph identity, we establish a lower bound for the convergence ra- 
dius which happens to be identical to the lower bound of the convergence radius of the virial 
series in the grand canonical ensemble established by Lebowitz and Penrose in 1964. We 
also show that the (Helmholtz) free energy can be written as a series in power of the density 
whose k order coefficient coincides, modulo terms o(N) /N, with the fc-order virial coefficient 
f*"*". ■ divided by fc + 1, according to its expression in terms of the m-order (with m < k + 1) simply 

connected cluster integrals first given by Mayer in 1942. We finally give an upper bound 
for the fc-order virial coefficient which slightly improves, at high temperatures, the bound 
obtained by Lebowitz and Penrose. 

o 

1 Introduction 

The rigorous approach to systems constituted by a large number of continuous classical interacting 
^ | particles has been a deeply investigated subject during the last decades. In particular, the study 

of the low density phase, where the behavior of the system should be near to that of an ideal gas, 
has yielded some of the most impressive results in mathematical physics. It is remarkable to note 
that most of these results have been obtained between 1962 and 1968. 

As far as a system (gas) of continuous particles is concerned, the Mayer series of the pressure 
(the pressure in powers of the fugacity) and the virial series of the pressure (the pressure in powers 
of the density) were known since the years around 1940. In particular, J. E. Mayer first gave the 
explicit expressions of the ra-order Mayer series coefficient in term of a sum over connected graphs 
between n vertices of cluster integrals, and of the n-order virial series coefficient in term a sum 
over two-connected graphs between n + 1 vertices of irreducible cluster integrals (see. e.g. [22] 
and references therein). However the question regarding the convergence of these series remained 
unanswered during the following two decades, since a "direct" upper bound of the type (Const.)" 
on these n order coefficients, which would have guaranteed analyticity of these series, was generally 
considered rather prohibitive, due to the fact that the number of connected (or even two-connected) 
graph between n vertices is too large (i.e. order C n with C > 1). 



1 



The rigorous analysis of the Mayer series and the virial series of the pressure of the system of 
particles started to produce results with the work by Groeneveld [18] in 1962, who first gave a bound 
of the type (Const.)" - for the n-order Mayer coefficient, but under the assumption that the pair 
potential is non-negative. One year later, Penrose [2S1 and independently Ruelle [3S1 ES] proved 
that the Mayer series of a system of continuous particles interacting via a stable and tempered pair 
potential (see ahead for the definition) is an analytic function for small values of the fugacity, as 
well providing a lower bound for the convergence radius. One year later Lebowitz and Penrose [115] 
obtained a lower bound for the convergence radius of the virial series (the pressure in function of 
the density). These results were all obtained "indirectly", (i.e. not by trying to bound directly the 
expressions of the coefficients as sum over connected graphs exploiting some cancelations) via the 
so called Kirkwood-Salszburg Equations (KSE), iterative relations between correlations functions 
of the system, and their possible use towards the control of the convergence of the Mayer series 
and virial series was glimpsed since the forties (see, e.g. |23|, I15j and reference therein). 

An alternative method to KSE was proposed in the same years by Penrose [30] who proved the 
convergence of the Mayer series of a system of particles interacting via a pair potential with a 
repulsive hard-core at short distance (but possibly attractive at large distance). To obtain this 
result he rewrote the sum over connected graphs of the n-order Mayer coefficient in terms of trees, 
by grouping together some terms, obtaining in this way the first example, as far as we know, of 
tree graph identity (TGI). It was only a decade later that Brydges and Federbush [5] were able 
to provide, for the second time, a proof of the analyticity of the Mayer series for the pressure of 
a continuous gas by directly bounding its n-term coefficients via a new type of TGI. Later, this 
direct approach based on TGI has been further developed and systematized by several authors 
(see e.g. [6J, [20], [5], [I], [35], [M], PU, 02]). Method based on KSE and TGI are part of the so 
called Cluster Expansion (CE) method. TGI are nowadays much more popular than the old KSE 
tools, mainly because of their flexibility and adaptability, but it is worths to remark that, as far 
as continuous particle systems interacting via a stable and tempered pair potential are concerned, 
the bounds given in [351 133 EFJ EH1 12] obtained via KSE have never been beaten. 

An other very popular tool in the framework of CE methods, which we need to mention here, is 
the so called abstract polymer gas (APG). The abstract polymer gas is basically a gas of subsets 
of some large set called polymers which posses a fugacity and interact via a hard core (non over- 
lapping) potential. Such model is in fact an extremely general tool for investigating analyticity of 
thermodynamic functions of virtually any kind of lattice system and its study has also a very long 
history which remounts to the sixties. Indeed the Polymer gas, as a gas of subsets of the cubic 
lattice "L d , was originally proposed in the seminal papers by Gallavotti and Miracle Sole [12] (as a 
tool for the study of the Ising model at low temperature) and Gruber and Kunz [13] (as a model 
in its own right), where the analyticity of the pressure of such gas at low densities is proved via 
KSE methods. The same model was then studied also using TGI methods e.g. by Seiler [UJ and 
Cammarota [7J. In 1986 Kotecky and Preiss [16] proposed a downright "abstract" polymer gas 
(polymers needed not to be subsets of an underlying set) and gave a proof of the convergence of 
the pressure not based on the usual KSE or TGI cluster expansion methods. They also provided 
criterion to estimate the convergence radius which improved those previously obtained via CE. 
The proof of the Kotecky-Preiss criterion was further simplified (and also slightly improved) by 
Dobrushin [5J[9] who reduced it to a simple inductive argument. The beautiful inductive approach 
for the polymer gas formulated by Dobrushin (which he called "no cluster expansion approach") 
was then generalized and popularized by the work of Sokal [ID] who also gave an extension of the 
APG convergence criterion for non-hard core repulsive pair interactions. Finally, the robustness of 
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Cluster Expansion has been recently revalued in [TO] where the Dobrushin Kotecky-Preiss criterion 
has been improved via TGI arguments [TO]. Moreover, always using TGI tools, further generaliza- 
tions of the APG convergence criterion for non-hard core, non-repulsive pair interactions has been 
recently given in [321 [33] , and [31] . 

Turning back to continuous particle systems, as mentioned above, all rigorous results about 
analyticity in the low density phase were obtained in the Grand Canonical Ensemble. Very recently 
Pulvirenti and Tsagkarogiannis [36] have obtained a proof of analyticity of the free energy of a 
system of continuous particles in the Canonical Ensemble. To get such a proof, they combined, 
within the cluster expansion methods, a standard Mayer expansion with the abstract polymer gas 
theory. In this work authors used the convergence criterion for the abstract polymer gas given in 
[2] and [25]. However, this is not the best criterion in the literature. In fact it is inferior not only 
to the recent Fernandez-Procacci criterion, but even to the Kotecky-Preiss criterion. Moreover, for 
technical reasons, in [36J the authors have used periodic boundary conditions (instead of the usual 
free boundary conditions), and consequently they had to assume some further condition on the pair 
potential beyond usual stability and temperness (see (2.3) in [36] and comments below). 

In this note we revisit the calculations in the canonical ensemble proposed by Pulvirenti and 
Tsagkarogiannis, but under the sole assumption that particles interact via a tempered and stable 
pair potential and are subjected to the usual free boundary conditions. Using the Fernandez- 
Procacci criterion to check the convergence of the polymer expansion, we show the analyticity of 
the free energy at low densities in the canonical ensemble and establish a lower bound for the 
convergence radius which improves the bound given in [36] and is identical to the lower bound of 
the convergence radius of the virial series in the grand canonical ensemble established by Lebowitz 
and Penrose [TO] in 1964. We also show that the ( Helmholtz) free energy can be written as a series 
in power of the density, whose k order coefficient coincides, up to terms o(N)/N, with the k order 
virial coefficient divided by fe + 1, according to its expression in terms of (simply) connected cluster 
integrals originally given by Mayer in 1942 (formula (49) in [21]). We finally give an upper bound 
for the /c-order virial coefficient which slightly improves, at least at high temperatures, the bound 
obtained by Lebowitz and Penrose in |19| . 

2 Notations and results 

Throughout the paper, if S is a set, then IS"! denotes its cardinality. If n is an integer then we will 
denote shortly [n] = {1, 2, . . . , n}. 

2.1 Continuous particle system: Notations 

We consider a system of N of classical, identical particles enclosed in a cubic box A C R d with 
volume V (and hence at fixed density p = N/V). We suppose N large (typically N 10 23 ). We 
denote by Xi G M d the position vector of the i th particle and by \xi\ is modulus. We assume that 
there are no particles outside A (free boundary conditions) and that these N particles interact via 
a pair potential V(xi — xj), so that the configurational energy of the N particles in the positions 
(xi, . . . ,xn) 6 h N is given by 





l<i<j<N 



We make the following assumptions on the pair potential V(x). 
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A. Stability: there exists B > such that, for all N £ N and for all {x\, . . . , xjy) € R rfAr , 

V{xi - Xj ) > -BN (2.1) 

l<i<j<N 

B. Temperness: 

C(P) = [ le-^W - l|da; < +00 (2.2) 



The (configurational) partition function of this system in the canonical ensemble at fixed density 
p = N/V and fixed inverse temperature (3 £ R + is given by the following function 



*a(/M = ^//*i-/ a < 



P Ei<i<j<w V(xi—Xj) 



dx N e P^i<i<i<N 

The thermodynamics of the system can be derived from the partition function Z\({3,p). In partic- 
ular, the Helmholtz free energy per unit volume of the system is given by 

/0M= lim /a CM (2.3) 

A,iV— >cx) 

where A — > 00 means that the size of the cubic box A goes to infinity and 

/aGM = --^1iiZaGM (2.4) 

We recall that the limit (2.3) is known to exists if the pair potential V(x) satisfies stability and 
temperness (see e.g. [37]). 

We will also make use in what follows of some notations concerning the system in the Grand 
Canonical Ensemble. We first recall the expression of the Grand Canonical Partition function of 
the system enclosed in the volume A, at fixed inverse temperature (3 and fixed fugacity A. 



Z A (p,\) = J2^J A d*i-J A dx N e- p ^<o<n 



with the N = term being equal to 1. The finite volume pressure of the system P\(f3,\) in the 
Grand Canonical ensemble is given by (see e.g. [TU] or |37| ) 

PP A (P, A) = i log H A (/3, A) = bn(P, A) A™ (2.5) 

n>l 

where b\{j3, A) = 1 and, for n > 2, 

&„(/?, A) = llf dxi ...f d x n J2 II [e~ mXl ~ Xj) " l] (2-6) 
n ' A A 9£G„ {i,j}eB g 

where G n is the set of all connected graphs with vertex set [n] and if g 6 G n then its edge set is 
denoted by The r.h.s. of (2.5) is known as Mayer series and the term 6„(/3,A) is the n-order 
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Mayer coefficient (a.k.a. n-order connected cluster integral). We will need in the next sections an 
upper bound for \b n (f3,A)\. This term has been bounded several times in the literature using both 
KSE or TGI methods. However, as far as we know, the bound obtained by Penrose nearly fifty 
years ago (see formula (6.12) in [28J) has never been beaten. The bound for |6 n (/3, A)| in [28] is, 
uniformly in A and for all n > 2, as follows 

|&n03,A)| < (2.7) 

re! 

where 

C(/3) = f le-W*) - l\dx (2.8) 



An easy computation (see again [19], [37]) shows that the finite volume density p = p\(/3,\) of 
the system in the Grand Canonical Ensemble is given by 

p = Y,nb n (f3,A)\ n (2.9) 

n>l 

So that one can eliminate A in (2.5) and (2.9) to obtain the so-called Virial expansion of the 
Pressure, i.e. the pressure in power of the density p = /?a(/3, A), in the Grand canonical Ensemble 

pp A (p, x) = P -J2 rryiW A )p k+1 ( 2 - 10 ) 
fc>i 

where, as shown more than fifty years ago by Mayer (see e.g. [22] and reference therein) 

P k (f3,A) = ~^dx 1 ---jjx k+ i II [e-^-^-l] (2.11) 



sgg* +1 {i,j}eE, 



with G* k+1 being the set of two-connected graphs with vertex set [k + 1]. The term /3fc(/3, A) is also 
known in the literature as the irreducible cluster integral of order k. 

As remarked in the introduction, the Mayer series in the r. h. s. of (2.5) has been proved to 
converge absolutely, uniformly in A |28|, [29"1 [38| [39] , for any complex A inside the disk 

' A ' < < 2J2 > 

where B is the stability constant defined in (2.1) and C{f3) is the function defined in (2.8). Moreover 
Lebowitz and Penrose [19] showed that the virial series in the r.h.s. of (2.10) converges for all 
complex p = pa(/3, A), uniformly in A, inside the disk 



with 



/ ^ \(l + u)e- w - l]w , 4 , 

g(u) = max ^ >- L_ 2 .14 

0<to<l U 



It is important to stress one again that the r.h.s. of (2.12), obtained in 1963, and the r.h.s. of (2.13), 
obtained one year later, still remain, as far as we know, the best lower bounds for the convergence 
radius of the Mayer series of the pressure and the convergence radius of the virial series of the 
pressure, respectively, of a system of continuous particles. 
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2.2 Results: a Theorem in the Canonical Ensemble 

To state the results obtained in this paper (resumed in Theorem 1 below), let us introduce the 
following notations. Following [36], we put 

Z A (j3,p) = ^-Z A (p,p) 

where 

dX! f dx^ N V(x l -x J ) 



Za(P,p)= Ht-- -^e- p ^<^<" v ^~ x >> (2.15) 
JA v Ja v 

Let us define 

QaCM = ^logZ A GM (2-16) 
Then the finite-volume Helmholtz free energy f\(/3,p) defined in (2.4) can be written as 

/aOM = -~ ^hx(^\+QMp) (2.17) 



where ^ln-^j- is the Helmholtz free energy of an ideal gas, and -^Qa(Pj P) is the part of the 
Helmholtz free energy due to the presence of the interaction V(x). 

Theorem 1 Let Qa(P,p) be defined as in (2.16), then the following statements are true, 
i) It holds that 

C fc (/3,A) 



Q^p) = E^jTT pk+1 (2 ' 18) 

k>l 



where, for any fixed k, 



n=l {m 2 ,...,m fc+1 } i=2 

u>ii/i i/ie 6j(/3, A) 's are i/ie (simply connected) cluster integrals defined in (2.6). 
ii) Let 

, _ -F(e 2 ^) 

P/J e 2 / 3B C(/3) 



(2.20) 



where C{(3) is the function defined in (2.8) and 



, . hl + ul - e" a 

Jm = max ) (| 2.21 

v ; a>o e a [l + u(l - e- a )\ 

Then the series in the r.h.s. of (2.18) converges absolutely, uniformly in A, in the (com- 
plex) disk \p\ < p* 3 . 



6 



Hi) As soon as the density p is such that p < p% the factors £fe(/3, A) in r.h.s. of (2.18) admit 
the bound, uniformly in A 



!<£*(/?, A) | < 



fc + 1 



+ e z 3 — lie z 3 



e 2/3B(fc-l) (fe + 1) [c , (/3)] fe (2 _ 22) 



loii/j being i/ie unique value of a £ (0, oo) suc/i i/iai i/ie function in the r.h.s. of (2.21) 
reaches its minimum value (i.e. reaches the value p%). 

Remark 1. The theorem above immediately implies that the infinite volume free energy f(/3,p) 
defined in (2.3) is also analytic in p in the same disk \p\ < p%. 

Remark 2. It is not difficult to check that J-(e 2B ^) is an increasing function of (3 with J-{1) ~ 
0, 1448 and lim^oo F{e 2B ") = e _1 . Moreover is a decreasing function of /3 with a^ =0 = 0, 426... 
and linig^oo a% = 0. We can thus compare this result with the best lower bound for the convergence 
radius of the virial series in the grand canonical ensemble (the pressure as a function of p). Such 
bound was obtained by Lebowitz and Penrose in 1964. They found that the virial series is analytic 
in the open disk \p\ < 1Z where (see formula (3.9) of [19J) 



with 



g{e 2pB ) 

U = e^ B C(P) (2 ' 23) 



/ ^ [(1 + u)e' w - l]w 
g(u) = max ± '- s — 2.24 

0<«J<1 u 



It has been for us quite surprising to realize that g{e 2 ^ B ) = J-(e 2 ^ B ), so that the lower bound 
(2.20) of the convergence radius of the Helmholtz free energy as a function of the density in the 
canonical ensemble and the lower bound (2.23) of the convergence radius of the virial series in the 
grand canonical ensemble given by Lebowitz and Penrose in [19] are in fact identical. 

Indeed, first note that the function G(w) = [( 1+M ) e ^ — zll^ inside the max in r.h.s. of (2.24) is 
positive only in the interval w € (0,ln(l + u)) and G(0) = G(ln(l + u)) = 0. So we can rewrite 

/ \ \(l + u)e- w -l]w i \ 

g(u) = max ^ '- s — (2.25) 

0<wj<1ii(1+m) U 



w \ 



On the other hand, via the change of variables w = ln[l + u{l — e a )] so that e a = u/(u + 1 — e 
the r.h.s. of (2.21) can be written 

ln[l + u(l -e _a )l x w(u+l-e w ) 

Fin) = max ) = max — - = g(u) 2.26 

v ; a>o e a [l + u(l - e- a )]' o<w<ln(i+«) e w u yx ' v ' 

Remark 3. Formula (2.19) is also quite remarkable since it immediately implies that £fc(/3, A) 

coincides, modulo terms of order ^9--, with the k order virial coefficient /3^(/3,A). Indeed, r.h.s. of 
(2.19) is, as was first shown by Mayer in 1942 [21] , the representation of the two-connected cluster 
integral /3jt(/3, A), defined in (2.11), in terms of the simply connected cluster integrals foj(/?,A) 
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(i = l,...,k + 1), defined in (2.6), (see [S] formula 49, see also formula (29) p. 319 of [SO]). 
Therefore we have, accordingly with [36J, that 



1 + 



o(N) - 
N . 



(2.27) 



Remark 4. We hnally compare our bound for the virial coefficients (2.22) with that obtained by 
Lebowitz and Penrose. Their bound, as stated in formula (3.13) of [19] is 



W/3,A) < 



0.28952 



(2.28) 



On the other hand bound (2.22) behaves asymptotically, taking for a*^ its largest (and hence worst) 
value a*p =Q = 0, 426, as 



Const. 



0.24026 



So our bound is asymptotically better than (2.28) for /3 small (i.e. high temperatures) and worst 
for (3 large (i.e. low temperatures). 

3 Proof of Theorem 1 



Following [36], we can easily rewrite Z\(f3,p) defined in (2.15) as a partition function of a hard 
core polymer gas. Indeed, using the Mayer trick we can rewrite the factor e~' 3 ^-' 1 ^ i< J^ JV v ( x <-~ x i) as 



'Ei< 



l<i<j<N 



g r- ^l<i<j<N " *3) = | | \e~' Jr _ I _|_ 1 j 



{R lt ...,Rs}eTT N 

where irjy = set of all partitions of [N] = {1,2,..., N}, and 



i 



if \R\ = 1 



E EI \e~P v ^- x ^ - 1] if \R\ > 2 

I gtG R {i,j}£E g 



with Gr being the set of connected graphs with vertex set R and, given g 6 2£ s denotes the 
set of edges of g. Now define, for any R C [N] such that \R\ > 2, 



<iHi = / A -/ A n^) 

JA i£R V 



(3.1) 



Note that depends only of the cardinality of the polymer R (the variables {xj}ie_R are mute 
variables). Observe also that 

b n (P,A)n\ 



yn-l 



(3.2) 
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where, for any n > 2, b n (/3, A) is the n order coefficient of the Mayer series of the pressure defined 
in (2.6). It is now easy to check that 



Za(P,P) = S[iv] = E E 



C]-Rl| • • • C\Rn\ 



(3.3) 



n>0 {«!,...,«„}: fliC[JV] 
|«j|>2, H in «j=0 



with the n = term giving the factor 1. So the partition function of a continuous gas in the 
canonical ensemble is equal to the hard core polymer gas partition function Enyi of a polymer gas 
in which the polymers are non overlapping subsets of the set [N] with cardinality greater than one, 
and a polymer R has activity given by (3.1). Note that the activity (mi can be negative since, 
by (3.2), its signal is the same of the Mayer coefficient &iffl(/3, A). 

It is well known (see e.g. [7], |35| ) that the logarithm of the hard core polymer gas partition 
function ^[n] can be- written as 



OO j 

lo s%] = E^ E 



n=l (Ri,..., J R„)G[7V] n 



(i?l, . . . ,Rn) (\ Rl \ ...C|i? n | 



(3.4) 



with 



( 1 



\Rl,...,Rn) = < 



96G„ 
^ 9 CG(il 1 ,...,il„) 



if n = 1 
-1)1^1 ifn>2 



(3.5) 



where G(R\, . . . , R n ) is the graph with vertex set [n] and edge set £-c7(iii,...,i?„) = {{*>j} ^ M : 
n ^ 0}, so Ylig&G n ls ^ ne sum over all connected graphs with vertex set [n] and U(R,R') 
which are also subgraphs of G(R\, . . . , R n ). In conclusion we can write the function Q\((3,p) 
defined (2.16) as 



OO 

3a08, p) = y E ^ E ■ ■ ■ > ^) <m ■ ■ ■ 



y — n! 

n=l (iii,...,i? n )e[JV] f - 



(3.6) 



3.1 Convergence criterion: proof of Theorem 1, part ii 

Once identity (3.6) has been established, or, in other words, once recognized that Qa(P,p) is the 
logarithm, divided by V of the partition function of a subset polymer gas (according with the 
terminology used in [3]), we are in the position to prove part ii of Theorem 1. By the Fernandez- 
Procacci criterion [10], we have that logEnyi defined in (3.4) (and hence Qa(P,p)) can be written 
as an absolutely convergent series for all complex activities C|_r| as soon as 



sup Yl K|/?iK |i?l < 

l ^l N ] RC[N]: ieR 
\R\>2 



(3.7) 



A few line proof of this statement can also be found in [3] (see there Theorem 2.4). Now, since the 
sum the l.h.s. of (3.7) does not depend on i £ [N] and using also the fact that the activity C|_r| 
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depends only on the cardinality of the polymer R, we can rewrite the condition (3.7) as 



A' 



e am C p m <e a -l (3.8) 

m=2 

where 

C p m =\Cm\( N ~]) (3.9) 
\m — 1/ 

An upper bound for the activity of a polymer R is now the key ingredient to implement the 
convergence criterion (3.7). Such a bound, recalling (3.2), follows immediately from the Penrose 
upper bound on the n-order Mayer coefficient given in (2.7). So we have 

I Cm I < V- m+1 e 2 ^ m - 2) m m - 2 [C{l3)] m - 1 (3.10) 

and, recalling that p = y, we get that C m admits the following estimate 

m-2 

C m < p — l e 2^(— 2)^L_^ [C(/3)r -l (3 . U) 

[m — lj! 

where we have bound {m-\)N~ m+1 < l/(m — 1)!. 
Hence, the convergence condition (3.8) is true if 

N m—2 

Y [pe-Ci^e 2 ^]" 1 - 1 ^ < e 2 ? B (l - er a ) (3.12) 

^-^ (m — 1)! 

m=2 v ; 



i.e. if 

n-1 



n 1 



n! 

n=l 



< l + e 2l3B (l -e- a ) (3.13) 



where k = l/(pe 2 ^ B C(/3)). Let now 



K* = min inf i k : Y ^ - < 1 + e WB (l 
a>0 [ ^ n\ k 



As shown in [3] (see also |llj and [2]) K* can be written explicitly as 

. e a [l+e 2 ^(l-e-°)] 
™>o ln[l + e 2 ^ s (l-e- a )] 

So (3.8), and hence (3.7), hold for all complex p as soon as 

ii^ * 
\P\ < Pp 

where 



and 



□ 



ln[l + u(l-e-°)] 
-Alii) = max — (3.15) 

v ' a>o e a [l + u(l - e- a )\ ' 
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3.2 Free energy in powers of the density: proof of Theorem 1 part i 

To prove formula (2.18), we reorganize the expansion (3.4), i.e. 



oo 

gaO?, p) = F E ^ E ^ T ( fl i. • • • > ^) ciflii • • • fa] 

n=l " (i?i,...,B n )e[AT]' 1 

as a power series in the density p. We use first of all the Penrose identity [30j 113 EH ESI E] which 
states (we use the notation of |10j ) that 



\R 1 ,...,R n ) = (-l) n - 1 W, 



G(fl 1> ...,fl„) 
r£T n 

where -Pg(Ri,...,.r„) are the set of Penrose trees of the graph G(R\, . . . , R n ) with vertex set [n] and 
rooted in a fixed vertex of [n], e.g., with root in the vertex 1. Thus 

Q^p) = vT / ^—H E w G( H 1 ,..., fln) Ci* 1 |---C|K n i (3.16) 

n=l rST„ (iJi,...,fl n )e[AT] n 

We put 

Cs = P*"V. (3-17) 

where, recalling (3.2), 

M/3,A)s! 



Then 



Qa(/3,P) = "TtE ~\ E E E 1 reP G{Rl ,..., Rn) C\R 1 \---C\R n \ 



N ^ n\ 

71=1 r&T n »l.".»n («! fl„)6[JV]>-- 

s> 2 



n E " E Csi-..Cs„E E ^e^GCHi,..,^) 

n=l ' »i.."«n T gT„ (Hi,....fln)e[iV]« 



n-1 



s,->2 



|JJi|=Sj 



n-1 



£ E t^r- E »" ~ V„ • • • V,„ £ E w ow „ 

n=l ' »i.-«n r( =T n (H 1 ,...,fl n )e[jv]i 



n-1 



= ^EM)r- E pV+i-pV+iE E Wg^,...,^) - 

71=1 fcl,.-.fcn TGT n («!,..., ii»)S[JV] n 

fcj>l \Ri\=ki + l 

^E — ~\ E^ E Mfe 1+ i---Mfc„+i E E 1 TeP G( H 11 „., Bn) 

71=1 fc>n fej.-fen: kj>l r£T n (.%,.. .,R»)e[iV]« 

k 1 +...+k n =k \R i \=k i +l 
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p 00 (-1)™- 1 k 
= nY1 -} E^ E /*«■■■/**, E E ^eP G{Rl ,..., Rn 

n=\ k>n si,...s„: Si >2 T eT n (R 1 ,...,R n )e[N]" 

si + ...+sn=k+n | J?^ | =s^ 



sE/E^T E E 



;i e^~)i — /'"j - ;; ~ 1^ u -•'■< .-. - 

fc>l n=l si,... Sn ; Sj >2 T £T n (fli,...,fl„)e[JV]" 

3 1 + ...+s„=fc+n |Hj|=3j 



^Ef^E^ E flKw.t-Mjj&T. E w e „, „ 

k>l Tl=l s 1 ,...s n :s i >2 j=l rGT n (fli,...,-Rn)e[JV] n 

s 1 + ... + s n = fe + n 1^1 = 8^ 

In conclusion we have proved formula (2.18) in Theorem 1, i.e. we have that 

Q^,P) = E^rrP k+1 ( 3 - 19 ) 



fc>i 



where 



and 



(r fc (/3,A)=^(-l)"- 1 ^ n (fc) (3.20) 



n=l 



7 . -I ™ 



»i, ■■■,*»» i=i 

Si =2,3,..., 
si + ...+s n =fc+n 



with 



^i, • • • . *n) = ^TT E E W^....,*., ( 3 - 22 ) 
reT n (Hi,...,-R„)e[Af]« 

|flll=si-.-,|-Rnl=Sn 

We now have to show formula (2.19). We start making the following observation. For fixed k 
and n, the integers si, . . . , s n such that si + . . . + s n = k + n define uniquely a fc-tuple of integers 
{rri2, • • • , w-fc+i} such that mj G {0, 1, 2, . . .} and {#j € [n] : Sj = i} = with the property that 
Ef=2^ = " and ^^(i - l)™, = fc and so U]=i[h ,(P, A) Sj \] = U^iW ^T 1 - Hence we 
can write 

k 4- 1 

^ n(fc) = ^r £ ntw^r* E n*. ■■■,*») (3.23) 

{mj = {m 2 ,...,m fc+1 } j = 2 S l;;- ! » 
T>i>0, Y, k itl m i= n {#je[n]?s~=i}=m i 
E l fc + 2 1 (<-l)m i =fc 

Let us now calculate P(si, . . . , s„) under the conditions that s±, . . . s n is an n-tuple of integers such 
that si > 2 and si + . . . + s n = k + n and such that it defines the fc-tuple of integers {7772, . . . ,777fc +1 }. 
Recalling thus (3.22), fix a tree r € T n and consider the factor 

E Wo^,...,^) (3-24) 
(j?. 1 ,...,n„)e[JV]™ 

|=s 1 ...,|i? 7l |=s n 
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This factor is clearly a polynomial in N and, in view of the limit TV" — > oo, we will retain only the 
term of maximal degree in N. It is not difficult to see that once the tree t € T n has been fixed, 
the contribution of the higher order in iV comes from the sum over polymers Ri submitted to the 
following prescriptions. 

• For any fixed vertex i of the tree r different from the root with degree di and i\, . . . , i^-i 
children sum over R( 1 , . . . Ri d ._ 1 in such way that each one of the polymers R^, . . . Ri d ._ 1 
shares exactly one vertex with Ri and all these di — 1 vertices are distinct (so in particular 
we must have di — 1 < \Ri\, otherwise the contribution of this tree vanishes). This is 
because in such way the factors 



N 



Ril I 



N 



R 



1 



have the maximal power in N. For the root 1 of r, with degree d% and i\, . . . , children 
do analogously (but this time it must hold d\ < \Ri\). 

• In any fixed vertex i of r we can choose the remaining \Ri . \ — 1 vertices in each polymer R4. 
associated to the ij child of i (j = 1, . . . , Si) among N vertices; actually the \Ri j \ — 1 vertices 
in each Ri. should be chosen among < N vertices, where N — k — n < < N — 1, 
according to the constraints imposed by the Penrose condition r € -Pg(.Ri,...,r„)> but, for k 

fixed, we have that k + n < 2k, so that N(l - 77) < N {j < N(l - i). 

These prescription are exactly the same conditions used to calculate the r.h.s. of equation 3.20 
in Lemma 3.4 of [llj. Proceeding thus analogously to the computation explained in Lemma 3.4 of 
[TT] we have, for any r 6 T n 



■10 T 



i=2 



N 
Si - 1 

i-i 



Si 

di - 1 



(di - 1)! 



1 + 



«i" \ d J ''fi^-l)! \di 



(di - 1)! 



1 + 



o(N) 
N 

o(N) 



^di<si,di<Si+l 



N 



di<si,di<Si+l 



Si! 



si! di\(s\ — d\ 



^■U 



1 



i=2 



( s i 



(di 



N 



A'+l 



1)! (di - l)\( Si - di + 1)! 

o(N) 



IV. 



1 + 



o(N) 
N 



di<s\,di<Si+l 



(si -di)\ Al ( Si 



n 



di + l)! 



1 + 



N 



II 



d\<si ,di<Si+l 



where di is the degree of vertices i in r and we recall that for any tree r £ T n we have that 
d\ + . . . + d n = 2n — 2. So, summing over all trees r 6 T n and recalling Cayley formula, we have 
that 



. . . ,s n ) 



E 



(n-2)! 



nr=i(*-i)!(si-^i)!H( S i-*+i)!L 



n 



1+ 



o(N) - 
N . 
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(n - 2)! 

(n - 2)! 



E 



(»!-!)! 



n 



Si! 



d 1 +... + d n =2n-2 
l<d 1 <s 1 , l<di<Sj + l 



(di - l)!(si - di)! (di - l)\( 8i -d t + 1)1 



1 + 



N 



nr=i(^-i)' 



E 



fr- 1 * 1 TT. 
- z, -in 11 



S,;! 



!!,...!„: !j>0 
il + ... + in=n-2 
0<Ji<si-l, 0<!j<s, 



iiKsi-Zi-lJIllliKsi-IOlL 



("-2)! 

_ nr=i(^-i)! 

In conclusion we have obtained 

(n-2)! 



E 



il,...i„: (;>0 
( 1 + ...+i„=n-2 
0<Ji<si-l, 0<!j< Si 



si-1 



n 

i=2 



1 + 



N 



1 + 



o(iV) 
N 



V(si,...,s n ) = 



nS[(i-i)ir 



E 



Il,...I„: i;>0 
il+... + !n=n-2 
0<i 1 <s 1 -l ; CKi^s; 



si-1 



n 

i=2 



1 + 



o(7V) 
N 



We now show the identity 



E 



*!,...(„: ij>0 
Jl + ...+l n =»-2 
0<i 1 <s 1 -l, (Kl^Sj 



n 

i=2 



k — 1 + n 
n-2 



(3.25) 



(3.26) 



Indeed, put t\ = s± — 1 and tj = Sj for all z = 2, . . . , n. So we need to prove that for any n-tuple 
ti,...,t n such that Y^=i U = n + k — 1, ii > 1 and tj > 2 for i > 2. 



e n ; 

il,...i„: !;>() 1 = 1 V 
il + ... + !n="-2 
0<i,<t; 



— 1 + n 
n-2 



(3.27) 



To prove this identity, suppose to have a set V with n + k — 1 objects and let V = V\ U V2 . . . U V n 
with V\,...,V n being n disjoint sets with cardinality \V\\ = t\,...,\V n \ = t n . Let K n ^ be the 
number of ways to pick up n-2 objects from the n + k — 1 objects of V. Of course K n ^ = ( fc ri^ n ) 
and this is the r.h.s. of (3.27). On the other hand, since the sets V±, . . . , V n form a partition of 
the set V we can also compute K n \. by choosing, for each i <G [n], b L objects from V; L which is 
done in ways, and then summing over all possible n-tuple fa, . . . , l n under the constraint that 
fa + . . . + l n = n — 2, getting in such way that K n ^ is also equal to the l.h.s. of (3.27). This prove 
(3.27) and hence (3.26). Putting this into (3.25) we have 



V(si, ... ,s n ) 



(n-2)! 



nSK*-!)!]™ 4 



fc-l + n\ o(N). 
n-2 J N ' 



Note that (3.28) implies that the factor P(si, . . . , s n ) is, at least modulo terms of order 
symmetric function of s±, . . . , s n , i.e. it depends only on n, k and numbers {m-2, . . . , nik+i}- Actually 



(3.28) 



N ' 
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one can easily realize that V(s%, . . . , s n ) is globally symmetric. Indeed, from definition (3.22) it 
immediately follows that V(si, . . . , s n ) is invariant under permutation of s%,. . . ,s n if the root of 
the Penrose trees is chosen to be 1. Moreover, by the construction of the Penrose identity, the 
r.h.s. of (3.22) do not depend on the choice of the root, so that in the end V{s\, . . . , s n ) is actually 
invariant under permutation of si, . . . , s n . 
We can now plug (3.28) into (3.23) and we obtain 

A; -1-1 

w n (k) = tt±\ 1+ °M.} y TT[W.A)iir J n ~ 2)! — ( k ~ l + n ) y i 

m i>°. £j= 2 m i= n {#j6W'»j=«}=»>i 
E^i 2 1 (i-l)m i =fe 

So, since 



we obtain 

w , w . (t+1)(n _ 2)![l+ ^)] E fimmi 

{ m i} = { m 2.---,™fc+l} i=2 

(*-l)"»j=fc 



£ i = 

Si >2 
{#jS[n]:s •=i}=m i 



(fc-l + n)! o(AQ v ffi fofoA)*] 

fc! iV ^ 11 mi! 

{ m i} = {™ 1 2v>™fc+i} 1=2 
Ej=2 ( 4 - 1 ) m i=' c 



and so 



WlA)=tl+ ^ 1E( - ir .(*4uoi E nJMMiir (3 , 9) 



71=1 {m i } = {m 2 ,...,m fc + 1 } 1=2 

Ej=2 f 4 - 1 )™^* 

which concludes the proof of part i) of Theorem 1. □ 

As previously remarked, r.h.s. of (3.29) is, up to terms of order ^zp-, exactly the expression 
given in formula (49) of [21] (see also (29) p. 319 of |26j). We can thus conclude that A) is, 

up to terms of order -^p, the very same k order virial coefficient as it is defined in formula (13.25) 
pag. 287 of |22j . So can also be written in terms of a sum over two-connected graphs between 

k + 1 vertices (see e.g. (13.25) in [22J). Namely, <£ fc (/3, A) = [1 + ^p-]f3 k (/3, A) where /3 k (b, A) is the 
virial coefficient defined in (2.11). 
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3.3 Bound for the free energy: proof of Theorem 1, part Hi 

Let us define the positive term series 



1 °° 1 

Ma(/M = ^£^ E E *TeP G(Rl ,..., Rn) \C\R 1 \\...\(\R n \\ (3-30) 

n=l ' r£T n (R 1 ,...,R n )e[N}" 

Then, clearly 

\Qa(P,A)\ < \Q\a(P,p) 



Now 



IQ|a(/3,p) = ^E^T E E E 1 reP G(Rl _ Rn) \C\R 1 \\---\C\R n \\ = 

n=l ' r€T„ n.^n (fl 1 ,...,.R ?l )e[.(V]™ 
S ^ 2 |Hil=»i 

= ^E^E^ E E E flri )Ki-ki = 

n=l fc>n si.-sn r gT n (fl 1 ,...,fl„)e[JV]" 

S1 +...s„ = fc+n |flj|=Sj 

= ^E^ fc E^ E E E w G(fll ,..., fl jM S1 i ■■>*„! = 

fc>l n=l ' »l.-»n T gT„ («i,...,il„)e[ivp 

sj + ...+% = Hn 



-Ef^E^ E nnw.*)M^E E w w „ „„, 

k>l n= ' • ' .. .. 

In conclusion we get 



k>l n=l »l.-*n i=l ' t£T„ (fl!,...,fl„)e[JV]« 



\Q\m P ) = p-^^ 1 P k + 1 (3.3D 

where 

fc 

\£ k \(P,A)=J2\W\n(k) (3.32) 

n=l 

with 

jU _|_ i n 

\W\ n (k) = -L- E nOMAA)Mn*i,.-.,Sn) (3-33) 

s 1 + ... + s n =fc + n 

Of course we have that 

|£ fe (/3,A)| < |Cfc|(/3,A) (3.34) 

We now obtain an upper bound for |Q|a(/3, p) as soon as p < p*^. We start by writing, recalling 
(3.30) 

1 °° 1 

iqu(m = f E^E^ ( 3 - 35 ) 

n=l " r€T n 
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where 



IC|Ri|l • • • \C\Rn\ 

(Rl,...,Rn)E[N] n 



We first observe that the n = 1 term in the sum in the r.h.s. of (3.35) can be written as 

777jE n - = 77 E IC|Ji|[ = ^Z}lC.| E 1 = ^E( J(^i = 



HC[JV] 
|fl|>2 



s>2 RC[N] 
\R\=a 



s>2 



1 N (N - V 

s>2 v 



1)101 =pE^ 

7 s>2 

Now following for a fixed n > 2 and r £ T„, we can bound 

IC|fli|l • • • \C\R n 



(3.36) 



Kl < E lr 

(fli,...,JS re )6[JV] 



CP G(R 1 ,...,R n ) 



where P, 



G(R u ...,R n ) 



are the weakly Penrose trees which are subsets of G(Ri, . . . , R n )- Using once 



again lemma 3.4 in [TT] we have, for any r E ^g(Ri r ) an( ^ an y n — 2, 



w< e w*iiC l A l ) ,K, n 

RlC[iV] \ J- / j_ 2 



«lC[iV] 
|flll>dlA2 



sup 

l«il>(di-l}A3 



l (di - mc ]Ri \ 



= w(di,. ..,dn) 

where di is the degree of vertex i of r. Hence, by Cayley formula and definition (3.9) 

(n-2)! 



1 x - 1 v - (n 



(* - 1)! 



«>(di,. • • ,d n ) 



d 1 + ... + d rl =2n — 2 



n(n — 1) 



e e K*i('*')*n 

,..-,dn BiClW v 17 i=2 



dl,...,d„ «iC[JV] 

d 1 + ... + d n =2n-2 |iJ 1 |>d 1 A2 



sup 

ie[iV] s,^,^ 

l-Kj|>(dj-i)Aa 



l-Ril 

d,- - 1 



\Cr, 



N 



e e #C:)*n 



n(n — 1) s\ \d\, 

y ' d x ,...,dn s 1 >d 1 f\2 V y 1=2 

d 1 + ... + d n =2n-2 



E 



s»>(d 1 -l)A2 



di-1 



CI 



we now use the trick first used in |35j (see there section 3), so that, multiplying and dividing by 
a n_1 (with a > 0), we get, for any n > 2, 



ni 



E i n -i 



< 



rGT„ 



n(n — 1) 



E E 

d n Sl>d\A2 
d 1 + ... + d n =2n-2 



_»i Mi 
si Vdi 



(iia dl J] 



i=2 



E 

S< >(di-1)A2 



di-1 



17 



< 



E E(2)*«*II 



E 

di-l=0 



Si \ — 1 



AT ^YP " 

— Y ^C§. ■■■CP 8ia(l+a) Sl - x T\[l + a\ Si 



N 



s>2 



n-1 



n(n-l) ^ 51 v 

V ; si>2 

Now, choosing a = e a P — 1, we get, by the convergence criterion (3.8), that for p < p*p 

i^(l + a) s Cf <1 



s>2 



So we get, for p < p*p 



1 EiEKi^^- 1 )E^ (sl " 1) 



F ^ n! 

n>2 t£T„ 



si>2 



and hence, (3.37) together with (3.36) yield the bound 



s>2 



Therefore, recalling bound (3.11) for C§ , we get 

\Q\a(P, P)<Y1 fxy [ 1 + ( k + l)(e°* " l)e' 



fc>l 



e 2^B(fe-l) ( fc + *)* 1 



By comparing (3.38) with (3.31) we immediately get 

1 



IW,A) < 



fc + 1 



+ (e a /3 - l)e 



alk 



(3.37) 



(3.38) 



and, by (3.34), the bound (2.22) follows. □ 

Acknowledgments 

Aldo Procacci has been partially supported by the Brazilian agencies Conselho Nacional de De- 
senvolvimento Cientffico e Tecnologico (CNPq) and Fundagao de Amparo a Pesquisa do estado de 
Minas Gerais (FAPEMIG - Programa de Pesquisador Mineiro). 



18 



References 

[1] A. Adbesselam and V. Rivasseau (1995): Tree forests and jungles: a botanical garden for 
cluster expansions in Constructive physics, Proceedings, Palaiseau, France 1994, Lecture 
notes in physics n. 446. 

[2] A. Bovier and M. Zahradnk (2000): A simple inductive approach to the problem of conver- 
gence of cluster expansions of polymer models. J. Statist. Phys. 100, 765-78. 

[3] R. Bissacot, R. Fernandez and A. Procacci (2010): On the Convergence of Cluster Expan- 
sions for Polymer Gases, J. Statist. Phys., 139, 598-617. 

[4] C. Borgs (2006): Absence of zeros for the chromatic polynomial on bounded degree graphs. 
Combin. Probab. Comput. 15, 63-74. 

[5] D. Brydges and P. Federbush (1978): A new form of the Mayer expansion in classical 
statistical mechanics . J. Math Phys. 19, 2064 (4 pages). 

[6] D. C. Brydges (1984): A short cluster in cluster expansions. In Critical Phenomena, Random 
Systems, Gauge Theories, Osterwalder, K. and Stora, R. (eds.), Elsevier, 129-83. 

[7] C. Cammarota (1982): Decay of correlations for infinite range interactions in unbounded 
spin systems. Comm. Math. Phys., 85, 517-28. 

[8] R. L. Dobrushin (1996): Estimates of semiinvariants for the Ising model at low tempera- 
tures. Topics in Statistics and Theoretical Physics, Amer. Math. Soc. Transl., 177, 59-81. 

[9] R. L. Dobrushin (1996): Perturbation methods of the theory of Gibbsian fields. In Ecole 
d'Ete de Probabilites de Saint-Flour XXIV - 1994, Springer-Verlag (Lecture Notes in Math- 
ematics 1648), Berlin-Heidelberg-New York, 1-66. 

[10] R. Fernandez and A. Procacci (2007): Cluster expansion for abstract polymer models. New 
bounds from an old approach, Commu. Math. Phys., 274, 123-140. 

[11] R. Fernandez and A. Procacci (2007): Regions Without Complex Zeros for Chromatic Poly- 
nomials on Graphs with Bounded Degree, Combin. Prob. Comp. 17, 225-238. 

[12] G. Gallavotti; S. Miracle-Sole (1968): Correlation functions for lattice systems, Commun. 
Math Phys. 7, 274-288. 

[13] C. Gruber and H. Kunz (1971): General properties of polymer systems. Comm. Math. 
Phys., 22, 133-61. 

[14] B. Jackson, A. Procacci and A. D. Sokal (2013): Complex zero-free regions at large \q\ for 
multivariate Tutte polynomials (alias Potts-model partition functions) with general complex 
edge weights, J. Combin. Theory, Series B, 103, 21-45. 

[15] J. G. Kirkwood (1946): The statistical mechanical theory of transport processes, J. Chem. 
Phys., 14, 180-201 

[16] R. Kotecky and D. Preiss (1986): Cluster expansion for abstract polymer models. Comm. 
Math. Phys., 103, 491-498. 



19 



[17] T. Kuna, Yu. G. Kondratiev, and J. L. Da Silva (1998): Marked Gibbs Measures via Cluster 
Expansion, Methods Funct. Anal. Topology 4, 50-81. 

[18] J. Groeneveld (1962): Two theorems on classical many-particle systems. Phys. Lett., 3, 
50-51. 

[19] J. L. Lebowitz and O. Penrose (1964): Convergence of Virial Expansions, J. Math. Phys. 
7, 841-847. 

[20] V. A. Malyshev (1980): Cluster expansions in lattice models of statistical physics and quan- 
tum theory of fields. Russian Mathematical Surveys, 35, 1-62. 

[21] J. E. Mayer (1942): Contribution to Statistical Mechanics, J. Chem. Phys., 10, 629-643. 

[22] J. E. Mayer and M. G. Mayer (1940): Statistical Mechanics, John Wiley & Sons, Inc. 
London: Chapman & Hall, Limited. 

[23] J. E. Mayer (1947): Integral equations between distribution functions of molecules, J. Chem. 
Phye., 15, 187-201. 

[24] S. Miracle-Sole (2000): On the convergence of cluster expansions. Physica A, 279, 244-249. 

[25] F. R. Nardi, E. Olivieri and M. Zahradnk (1999): On the Ising model with strongly 
anisotropic external field, J. Statist. Phys. 97, 87-144. 

[26] R. K. Pathria and P. D. Beale (2011): Statistical mechanics, Third edition, Elsevier, Ams- 
terdam. 

[27] Ch.-E. Pfister (1991): Large deviation and phase separation in the two-dimensional Ising 
model, Helv. Phys. Acta, 64, 953-1054. 

[28] O. Penrose (1963): Convergence of Fugacity Expansions for Fluids and Lattice Gases, Jour- 
nal of Mathematical Physics 4, 1312 (9 pages). 

[29] O. Penrose (1963): The Remainder in Mayer's Fugacity Series, J. Math. Phys. 4, 1488 (7 
pages). 

[30] O. Penrose (1967): Convergence of fugacity expansions for classical systems. In Statistical 
mechanics: foundations and applications, A. Bak (ed.), Benjamin, New York. 

[31] S. Poghosyan and D. Ueltschi (2009): Abstract cluster expansion with applications to sta- 
tistical mechanical systems, J. Math. Phys. 50, no. 5, 053509, (17 pp). 

[32] A. Procacci (2007): Abstract Polymer Models with General Pair Interactions, J. Stat Phys., 
129, 171-188. 

[33] A. Procacci (2009): Erratum and Addendum: "Abstract Polymer Models with General Pair 
Interactions", J. Stat. Phys., 135, 779-786. 

[34] A. Procacci, B. N. B. de Lima and B. Scoppola (1998): A Remark on High Temperature 
Polymer Expansion for Lattice Systems with Infinite Range Pair Interactions, Lett. Math. 
Phys., 45, 303-322. 



20 



[35] A. Procacci and B. Scoppola (1999): Polymer gas approach to N-body lattice systems. J. 
Statist. Phys. 96, 49-68. 

[36] E. Pulvirenti and D. Tsagkarogiannis (2012): Cluster Expansion in the Canonical Ensemble, 
Comm. math Phys. 316, Issue 2, pp 289-306. 

[37] D. Ruelle (1969): Statistical mechanics: Rigorous results. W. A. Benjamin, Inc., New York- 
Amsterdam. 

[38] D. Ruelle (1963): Correlation functions of classical gases, Ann. Phys., 5, 109-120. 

[39] D. Ruelle (1963): Cluster Property of the Correlation Functions of Classical Gases, Rev. 
Mod. Phys., 36, 580-584. 

[40] A. D. Sokal (2001): Bounds on the complex zeros of (di) chromatic polynomials and Potts- 
model partition functions, Combin. Probab. Comput. 10, 41-77. 

[41] E. Seiler (1982): Gauge Theories as a Problem of Constructive Quantum Field Theory and 
Statistical Mechanics, Lecture Notes in Physics 159, Springer- Verlag, Berlin-Heidelberg- 
New York. 

[42] D. Ueltschi (2004): Cluster expansions and correlation functions. Mosc. Math. J. 4, 511- 
522. 



21 



